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ABSTRACT 

The HIDEEP survey (Minchin et al. 2003) was done in an attempt to find objects 
, having low inferred neutral hydrogen column densities, yet they found a distribution 

J> ■ which was strongly peaked at 10 20 65 cm~ 2 . In an attempt to understand this distor- 

ts — ' bution and similar survey results, we model HI profiles of gas discs and use simple 

t**»« . simulations of objects having a wide range of HI properties in the presence of an ioniz- 

ing background. We find that inferred column density (Nfjj) values, which are found 
by averaging total HI masses over some disc area, do not vary strongly with central 
column density (N max ) for detectable objects, so that even a population having a wide 
range of N max values will give rise to a strongly peaked distribution of N^j values. We 
find that populations of objects, having a wide range of model parameters, give rise to 

■ inferred column density distributions around 10 20 6±0 - 3 cm~ 2 . However, populations of 

fairly massive objects having a wide range of central column densities work best in re- 
producing the HIDEEP data, and these populations are also consistent with observed 

■ Lyman limit absorber counts. It may be necessary to look two orders of magnitude 
^ \ fainter than HIDEEP limits to detect ionized objects having central column densities 

< 10 20 cm~ 2 , but the inferred column densities of already detected objects might be 
lower if their radii could be estimated more accurately. 



Key words: galaxies: structure, mass function, ISM: general, diffuse radiation, radio 
lines, intergalactic medium 



1 INTRODUCTION 

Understanding the properties of dwarf galaxies, large dif- 
fuse galaxies, and any clouds of similar mass is important in 
understanding the formation of galaxies. For example, Cold 
Dark Matter theory predicts the existence of a population 
of low-mass satellite galaxies (e.g. Moore et al. 1999; Klypin 
et al. 1999). Furthermore, studying the properties of such 
objects is important in understanding the nature of Lya 
absorbers and metal line absorbers such as weak Mgll sys- 
tems (Rigby, Charlton, & Churchill 2002). Any such objects 
which have yet been undetected may have a different range 
of averaged neutral hydrogen column densities from that of 
the known population of galaxies. 

Gas having a wide range of neutral column density 
(Nhi) values has been observed as Lya absorption at low 
redshifts (for example, Bahcall et al. 1996) where absorp- 
tion lines shortward of Lya emission in quasar spectra arise 
from lines of sight through intervening gas between us and 
the quasar, and Nhi ranges from < 10 12 cm -2 to ~ 10 21 
cm -2 . Larger amounts of gas with Nhi J; 10 19 cm -2 can 



also be observed more directly as 21 cm emission in the lo- 
cal universe. The strongest 'damped' Lya absorbers, with 
Nhi > 10 20 ' 3 cm -2 , are often found to arise in lines of 
sight through galaxies including several low surface bright- 
ness (LSB) and dwarf galaxies (Cohen 2001; Turnshek et 
al. 2000; Bowen, Tripp, & Jenkins 2001). Yet the somewhat 
weaker Lyman limit systems (Nhi > 10 17 ' 2 cm -2 ), the col- 
umn densities of which are more difficult to measure ac- 
curately, have long been thought to arise in lines of sight 
through luminous galaxies (Bergeron & Boisse 1991; Stei- 
del 1995). Some weaker Lya forest absorbers are thought to 
arise in small amounts of intergalactic gas (Dave et al. 1999), 
while some could arise in gas surrounding galaxies (e.g. Chen 
et al. 2001; Linder 1998; 2000). 

A recent HI survey (HIDEEP; Minchin 2001; Minchin et 
al. 2003) was capable of detecting objects with inferred neu- 
tral hydrogen column densities (N H i) as low as 4x 10 18 cm -2 
for galaxies having velocity width AV = 200 km s -1 , assum- 
ing that a galaxy with suitable properties fills the telescope 
beam. Yet they failed to find anything with N HI < 10 20 
cm -2 . Other HI surveys have also found that galaxies show 
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little variation in column densities averaged over some ra- 
dius (Zwaan et al. 1997), although the integration times may 
not be long enough to detect low column density galaxies in 
such surveys, as discussed by Minchin et al. (2003). These 
HI surveys are limited by flux, rather than column density, 
when detecting faint objects, and the column density of the 
detected objects is uncertain given that the sources are gen- 
erally unresolved. However there is a limit on column density 
in a survey such as HIDEEP in the sense that a resolved, 
low column density object could fill the beam, although such 
objects are not often seen. 

Rosenberg & Schneider (2003) found that their sample 
of Hl-selected galaxies obey a relationship between HI cross 
section and HI mass, which is equivalent to having fairly 
constant averaged column densities. They plot, in their first 
figure, the disc areas Aola, where Nhi > 2 x 10 20 cm~ 2 
and thus where damped Lya absorbers can arise, versus 
the HI mass (Mhi) for a sample of HI selected galaxies. 
Some scatter is seen in the log-log plot, yet they can eas- 
ily fit a line having a slope of about one. Thus they find 
\og(AoLA) = log(Affrj) — 6.82, which would imply that 
galaxies having a wide range of mass and HI sizes all have 
area-averaged column densities of around 8 x 10 20 cm -2 , 
where the displayed points are all within about 0.8 orders 
of magnitude from the fitted line. 

Similar correlations between HI size and HI mass have 
also been seen by Giovanelli & Haynes (1983) and Verhei- 
jen & Sancisi (2001), and a correlation between HI mass 
and optical sizes of galaxies has also been seen by Haynes & 
Giovanelli (1984). Other surveys, capable of detecting low 
HI mass objects at various sensitivities, including some di- 
rected toward detecting extragalactic high velocity clouds 
(HVCs) (Blitz et al. 1999; Charlton, Churchill & Rigby 2000; 
Davies et al. 2002) have been largely unsuccessful at finding 
objects with low HI masses (de Blok et al. 2002; Zwaan & 
Briggs 2000; Dahlem et al. 2001; Zwaan 2001; Verheijen et 
al. 2000). On the other hand, some very faint optical sources 
have been found to be rich in gas (Davies et al. 2001), and 
there is theoretically no reason to expect every HI cloud to 
be capable of forming large amounts of stars. Furthermore, 
small HVCs with peak Nhi ~ 6 x 10 18 cm -2 are being de- 
tected around our Galaxy (Hoffman, Salpeter, & Pocceschi 
2002) and around M31 (Thilker et al. 2004). 

One suggested explanation for the lack of low column 
density detections in the HIDEEP survey is that the gas is 
hidden in 'frozen discs' (Minchin et al. 2003) where the 21 
cm transition is not excited to a spin temperature above the 
cosmic background (Watson & Deguchi 1984). 

A second possible explanation for the lack of low col- 
umn density detections is that the gaseous discs become 
highly ionized at a disc radius not far beyond that where 
Nhi = 10 20 cm -2 , so that the average inferred value re- 
mains above 10 20 cm -2 . The ionization of outer galaxy discs 
by a background of Lyman continuum photons was sug- 
gested and modelled first by Bochkarev & Sunyaev (1977) 
and later by Maloney (1993), Dove & Shull (1994a), and 
Corbelli & Salpeter (1993) in order to explain the sudden 
truncations seen in carefully observed spiral galaxy discs. 
Since then ionized gas has been detected in Ha emission us- 
ing a Fabry-Perot 'staring technique' (Bland-Hawthorn et 
al. 1994) beyond the HI edges of several nearby galaxies 
(Bland-Hawthorn, Freeman & Quinn 1997; Bland-Hawthorn 



1998). The ionizing background has been measured most 
recently at low redshifts by Scott et al. (2002), who find 
J(912A) = 7.6t'tt x W 23 erg cm -2 s _1 Hz -1 sr" 1 . Gas in 
the ionized parts of outer galaxy discs is likely to give rise to 
at least some Lya absorption (Linder 1998; 2000), and some 
variations will arise in the column density value at which HI 
discs are truncated as a result of fluctuations in the ionizing 
background radiation (Linder et al. 2003). 

Ionized gas clouds cannot correctly be referred to as 
undetected 'HI clouds' (although current HI structures may 
have been ionized in the recent cosmological past). However 
structures containing ionized gas are interesting and rele- 
vant to the galaxy formation process. For example, ionized 
gas contains enough neutral atoms to give rise to all of the 
Lyman alpha absorbers (except for the damped ones), and is 
thus, in principle, detectable in deep HI observations. HVCs 
may also contain mostly ionized gas. It is unknown whether 
massive clouds exist far from luminous galaxies, although 
not all of the absorbers arise close to galaxies (Stocke et 
al. 1995). Ionized gas clouds may have small regions con- 
taining HI clouds if the gas is sufficiently clumpy. 

In this paper, we wish to understand the observed lower 
limits in averaged column densities, and to constrain the 
properties of any objects that could be going undetected in 
HI surveys as a result of photoionization. Section 2 discusses 
the modelling of HI discs and calculation of column densities 
from HI observations. Section 3 describes the method used to 
model galaxy and cloud HI profiles and simulate populations 
of objects having a wide range of properties in the presence 
of an ionizing background. The results of such simulations 
are discussed in Section 4. The value for the Hubble constant 
is assumed to be Ho = 80 km s _1 Mpc -1 . 



2 HI DISCS AND AVERAGED COLUMN 
DENSITIES 

Suppose all galaxies have exponential HI column density 
profiles (as modeled for example, by Swaters et al. 2002) 
with central column density N max and scale length h so that 
Nni{r) — N max exp(—r /h) at a radius r along the disc. If 
an observer maps the profile out to column density N m i n , or 
an equivalent radius of h \n(N max / N m i n ), the HI mass could 
then be found within an area A = irh 2 In 2 (N max / N m i„) . The 
HI mass within this radius, where Nhi > N m in, would be 



Mhi = 2irh m H N„ 



N„ 



N„ 



1 + ln 



N„ 



Nn 



(1) 



where van is the mass of a hydrogen atom. The averaged 
column density (Nhi) = Mhi /(thhA) thus becomes 



{Nhi) = 2N n 



In 2 



(2) 



which depends only on N max for a sample of galaxies 
mapped out to a constant N m i n , and not on the scale length 
h. However, a fairly constant (Nhi) does not imply a con- 
stant N max , as seen in Fig.0 where (Nhi) is plotted versus 
N max , as in equation (2), for several values of N m i n . A nar- 
row range of (Nhi) allows for a somewhat wider range of 
N max and only implies that N max values are not likely to 
be in the higher range shown on the plot, where the curves 
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Figure 1. Values of the average column density (Nhi) are plot- 
ted versus central column density N max , as in equation (2), for 
several values of limiting N m in . The curves are fairly flat for lower 
values of N max . Thus fairly constant observed column densities 
would allow for N max to have a wide range of moderate to low 
values. The formula (2) plotted here is valid only for purely ex- 
ponential discs, without considering ionization effects. 
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Figure 2. Values of inferred column density N^j are plotted 
versus central column density N m ax, as in Equation (3), again for 
exponential discs without yet considering ionization, for several 
values of limiting N m i„. Here the curves turn up for lower values 
of N m ax when a constant N m in is assumed, as most of the mass 
will be outside the assumed radii for these objects. Thus there 
is a minimum N^j that can be measured for a given N m in, and 
lower values can only be found if the radii are estimated 

differently for objects having low N max values, for example as 
shown in the line where log N m in = log N m ax — 0.5. 



become steeper. (We also know that damped Lya absorbers 
are not seen with Nhi <; 10 22 cm" 2 .) 

Some scatter is seen in the column density measure- 
ments from Minchin et al. (2003) and Rosenberg & Schnei- 
der (2003), and in central column densities found in studies 
of LSB galaxies such as de Blok, McGaugh, & van der Hulst 
(1996). Thus some variation is likely to occur for values of 
N max , though its distribution is unknown. 

Note that while the Rosenberg & Schneider (2003) 
data can be compared with the {Nhi) values shown above, 
Minchin et al. (2003) calculate column densities somewhat 
differently, so we will refer to these as 'inferred' column 
densities (N HI ) in order to distinguish from the averaged 
ones discussed above. In this case the complete HI mass 
of a galaxy is measured, assuming that the mass is entirely 
within the beam. Since they have not resolved their detected 
objects or obtained column density profiles, they use an es- 
timate of the HI radius in order to find an inferred column 
density. They assume that the HI radii are equivalent to 5 
times their observed effective optical radii, based upon the 
relationships from Salpeter & Hoffman (1996). Here we gen- 
erally assume that such a radius corresponds to that where 
Nhi = 10 20 cm" 2 , although Salpeter & Hoffman (1996) base 
their relationships on several studies which may have slightly 
varying limiting values. Since the N HI values are found dif- 
ferently from {Nhi) values, they will depend here on what 
fraction of the HI mass is contained within the estimated 
radius. The inferred column density for an exponential disc 
where the complete mass is accurately measured would be 

AT o _ Nmax / Q \ 

HI ~~ t 2/ at Tat V 1 ^' 

In {Nmax /Nmin) 

In this case N m i„ is simply the column density correspond- 
ing to the estimated galaxy radius. Like equation (2), this 
formula does not depend upon the disc scale length. 



Objects having Nmax ~ N m i„ will have more of their 
mass outside the radius at Nmin, which will cause N H i to 
become large as seen in Fig. |5] Thus there is a minimum 
N H i value that can be detected for a given constant N m in- 
Objects having Nmax ~ N m in are likely to have small radii 
however, and these radii may be difficult to measure in a 
manner consistent with that used for larger objects. Thus 
an observer might effectively be using a smaller N m in for 
a smaller N max when estimating any measurable radius for 
such an object. For example if \ogN m in = logjV m<M — 0.5 
then N H i = 0.745A r mal for an exponential disc, as shown 
by the dotted line in Fig. |5] 

Different galaxy radii are used in different studies, in- 
side of which the column densities are averaged or inferred. 
The data from Rosenberg & Schneider (2003) can be used 
to find averaged column densities for their sample out to a 
radius where Nhi = 2 x 10 20 cm -2 , but their plotted HI 
masses and disc areas do not tell us about the properties of 
galaxies whose central column densities are about this value 
or lower. Based on their fit to the points in their first figure 
with a slope of one, their characteristic average column den- 
sity would be (Nhi)* = 8.26 x 10 20 cm" 2 or 10 20 92 cm" 2 , 
which would correspond to N max — 10 21 ' 7 cm -2 accord- 
ing to equation (2). However, points scattered within about 
half an order of magnitude of this value would suggest that 
galaxies exist that have N max values well below 10 21 cm -2 . 
Many such objects are likely to have small HI radii above 
2 x 10 20 cm" 2 , and thus need to be detected in a deeper 
survey such as that of Minchin et al. (2003). 

Minchin et al. (2003) report a characteristic Nhi — 
10 20 ' 6 cm" 2 for the HIDEEP sample. If we assume the 
N max value found above, then Equation (3) would give 
N m i n = 10 21 ' 6 cm -2 , which would not allow for the de- 
tection of objects with very low, or even moderate, Nffi- 
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This N m i n value is probably too high to make sense or al- 
low for the detection of many objects. However the HIDEEP 
survey is likely to detect objects with lower typical Nmax val- 
ues than those of Rosenberg & Schneider (2003). Suppose 
we assume that the radii used by Minchin et al. (2003), who 
use five times the effective optical radii, are equivalent to 
those where Nmin ~ 10 20 cm -2 as is typical of the galaxies 
studied in Salpeter and Hoffman (1996). In this case, Equa- 
tion (3) would give two solutions: N ma x = 10 2 cm -2 and 
N m ax = 10 20 ' 3 cm -2 , as seen in Fig. [5] The larger value 
would be unusually high compared to what is seen in galax- 
ies with known HI profiles or damped absorbers, although 
the curve in Fig.|5|is not very steep for high N max values. On 
the other hand, objects having the lower N max value would 
likely have small radii so that the radii could be estimated 
less carefully than for larger objects. It is thus likely that 
objects with a wide range of N max values, including lower 
ones, are being detected by HIDEEP. 

A larger radius, where N m in < 10 20 cm -2 , would be 
needed to calculate N^j below 10 20 cm -2 . However the 
limiting parameter for detecting clouds and galaxies in HI 
surveys is flux, rather than column density. Thus a more 
interesting question, as opposed to understanding the accu- 
racy of galaxy radii, is whether objects having low values of 
N m ax and sizes or masses similar to those of known galaxies, 
could be detected in a survey having some limiting flux. We 
discuss this issue further using simulated galaxies in Sec. 4. 

Note that column density profiles may fall off more 
slowly than exponentials in the outer parts of galaxies (Hoff- 
man et al. 1993; see discussion in Linder et al. 2003). Fur- 
thermore, exponential profiles are not always well-behaved 
in the centres of galaxies where there may be stars or molec- 
ular gas instead of neutral hydrogen. (Thus N max becomes 
the extrapolated central column density assuming an expo- 
nential profile in the galaxy's centre.) Finally, column den- 
sity profiles are thought to fall off quickly at a few x 10 19 
cm -2 due to ionization, which was not yet considered in 
Figs.0and|51 The effects of using more realistic column den- 
sity profiles will be discussed further for simulated galaxies. 



3 HI PROFILES AND SIMULATIONS 

Simulations are done in order to determine the HI fluxes 
for possible populations of low column density objects, and 
to attempt to reproduce the distribution of inferred column 
densities seen in the HIDEEP survey. Samples of galaxies (or 
clouds, having unknown optical properties) are simulated at 
z = in order to produce the figures discussed in Section 
4. Gas in each galaxy is modelled as a slab structure in 
hydrostatic equilibrium, where the gas is confined by a com- 
bination of pressure and gravity as in Charlton, Salpeter, & 
Hogan (1993) and Charlton, Salpeter & Linder (1994). 

We wish to simulate objects having a wide range of 
properties, especially those which are low in mass or col- 
umn density which may be difficult to detect due to ioniza- 
tion. We assume that each galaxy has an exponential total 
(neutral plus ionized) column density profile to start. Fur- 
ther simulations vary the profile, for example using a power 
law fall-off beyond four HI disc scale lengths, as discussed 
in Linder et al. (2003). The central (total) column density 
N m ax is assumed to have values which are uniformly dis- 



tributed between 10 18 and 10 22 ' 2 cm -2 . The higher central 
column density limit is chosen as an upper limit of what is 
seen in detected galaxies and damped Lya absorbers, while 
the lower range is used in order to explore the possible exis- 
tence of gas clouds which are more difficult to detect, being 
somewhat below the current sensitivity of HI surveys. 

Disc scale lengths are chosen so that the simulated ob- 
jects obey a Schechter-type total gas mass function, which 
gives rise to a detectable HI mass function having a similar 
slope of -1.3 (Zwaan et al. 2003) or -1.5 (Rosenberg & Schnei- 
der 2002). Two main cases are illustrated repeatedly in the 
following section: In Case A, we assume that the central 
column densities of the objects are correlated with the total 
gas masses. Case A is motivated by HI observations of LSB 
galaxies which suggest that they have lower central column 
densities (de Blok et al. 1996) and the likely existence of nu- 
merous dwarf LSB galaxies, such as in Sabatini et al. (2003) . 
For each simulated galaxy a relationship is assumed for Case 
A where log N max = 21.7+1.0 log(A/ tot - 10.0) ±0.5, so that 
objects having the lowest simulated N m ax values will tend 
to also have the gas masses around 1O 6 M0. (We later vary 
this relationship, as the narrow range of scatter is chosen as 
an extreme example to start.) In this case we have a sub- 
stantial population of small clouds having low N ma x, some 
of which might resemble HVCs, although the detectable ob- 
jects will tend to have high N m ax- In Case B the disc scale 
lengths and the central column densities are uncorrelated. 
Thus N max is uniformly distributed and unrelated to the 
total gas mass, but the scale lengths tend to be larger for 
objects having lower N max . In this case we are simulating 
larger objects having low column densities, which might re- 
semble giant LSB galaxies (whose numbers are very uncer- 
tain) or extended structures that give rise to Lya absorption, 
as expected based upon double line of sight observations (for 
example, Dinshaw et al. 1998; Charlton, Churchill & Linder 
1995; Monier, Turnshek & Hazard 1999; Fang et al. 1996). 

Rotation velocities are found for each simulated galaxy 
using the relationships given in Salpeter & Hoffman (1996), 
where the Tully-Fisher relationship between the veloc- 
ity Vrot and the observable HI radius R is found to be 
Kot/80.51 kms" 1 = {R/12.3 kpc) 1 ' 38 . The value of R is 
assumed to typically correspond to a radius where the lim- 
iting column density is 10 20 cm -2 , but then this relationship 
does not give us information about the rotation velocities of 
massive objects having N m ax < 10 20 cm -2 . Since galax- 
ies having a wide range of properties are found to obey a 
baryonic Tully-Fisher relation (McGaugh et al. 2000), we 
extrapolate the HI Tully-Fisher relation above into a bary- 
onic version by assuming that R in the formula above is the 
radius that a galaxy of equivalent (neutral plus ionized hy- 
drogen) mass would have if it had N max = N ma x, where we 
assume N ma x = 10 21 ' 7 cm -2 as found in the previous sec- 
tion. When we attempt to simulate the smallest clouds, the 
value of R may be very small, so that the velocity disper- 
sion of the gas becomes more important. A minimum value 
of Vrot = 10 km s _1 is thus assumed. Note, however, that 
there is a selection effect against detecting objects having 
velocity widths <J 50 km s _1 in HI surveys, as discussed in 
Minchin et al. (2003) and Lang et al. (2003). 

The vertical ionization structure of the gas is modelled 
as in Linder (1998), which is similar to the model in Mal- 
oney (1993). Inside of some ionization radius R cr , the gas 
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is assumed to have a sandwich structure, where the inner 
shielded layer remains neutral and has a height (zi) de- 
termined by equation (6) in Linder (1998). The gas above 
height (zi) and beyond the ionization radius is assumed to 
be in ionization equilibrium. 

The frequency- and direction-averaged ionization rate 
£ is assumed at first to be 3.035 x 10 -14 s _1 , from the cal- 
culation of Dave et al. (1999) at z — based upon spectra 
from Haardt & Madau (1996). The lowest measurements at 
redshifts ~ tend to be consistent with this value as dis- 
cussed in Linder et al. (2003). However galaxies, in addition 
to quasars, may contribute to the ionizing background radi- 
ation (Giallongo, Fontana, & Madau 1997; Shull et al. 1999; 
Bianchi, Cristiani, & Kim 2001; Linder et al. 2003). The 
ionizing intensity may be stronger when close to a luminous 
galaxy or galaxy-rich environment as shown in Linder et 
al. (2003), although the gas-rich objects which we are sim- 
ulating here are not likely to arise in the most galaxy-rich 
environments. Thus simulations are also run using a larger 
frequency and direction-averaged ionizing intensity measure- 
ment for redshifts z < 1 of ( = 1.9 X 10" 13 s" 1 (Scott et 
al. 2002). The conversion between £ and a one-sided flux 
is assumed as in Tumlinson et al. (1999). The radius R cr , 
at which the disc becomes fully ionized, is found where the 
neutral gas height Zi becomes zero. 

For each galaxy, the neutral column density profile is 
calculated by integrating the neutral density uhi verti- 
cally through the disc, in increments of scale length h/10 
(or smaller when needed for a more accurate mass calcula- 
tion). The profiles have Nhi ~ Ntot in the regions inside of 
Rcr- Just inside the radius R cr , the column density falls off 
quickly, typically from Nhi ~3x 10 19 cm" 2 to Nhi ~ 10 17 
cm -2 , at which point the 'ionized region' of the disc is being 
mapped, and only a small fraction of the gas is neutral. The 
resulting profiles can then be integrated and averaged over 
suitable radii to be compared with HI observations. 

In order to calculate fluxes for the simulated objects, 
each object is assigned a random inclination and a random 
distance within a sphere around us having a radius of 108 
Mpc, the distance at which a 10 M© galaxy can be detected 
at a limiting peak HI flux of 18 mjy/beam, as in Minchin 
et al. (2003). The HI mass for each object, limited to what 
can be contained within a 15 arcmin beam centred on the 
galaxy, is binned into velocity channels of 15 km s _1 when 
finding a peak or integrated flux, which is comparable to 
what is done for the HIDEEP survey. 



4 RESULTS 

Averaged column densities are plotted, again versus cen- 
tral column density N m ax, for simulated galaxies, which are 
exposed to an ionizing background. For simulated galaxies 
N ma x is defined as the total (neutral plus ionized) assumed 
central column density, which is about equal to the observ- 
able neutral value for N max > a few xlO 19 cm -2 . For lower 
Nma X values, the neutral values could be as low as a few 
xlO 17 cm -2 , although it is difficult to determine the value 
accurately very close to the centres of the discs with the 
model used here. Gas with column densities <, 10 19 cm -2 
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are seen, for example as mini-HVCs (Hoffman et al. 2002). 
Purely exponential total column density profiles are as- 



Figure 3. Averaged neutral column densities are plotted ver- 
sus central neutral column densities (N max ) for simulated galax- 
ies having exponential column density profiles with h = 2 kpc 
which are exposed to an ionizing background. Column densities 
are found here by integrating HI masses out to N m in and dividing 
by the area within the radius of N m i n , where N m in ranges from 
2 X 10 20 cm -2 (top curve) to 10 18 cm -2 (bottom), which is indis- 
tinguishable from N m i n = 10 19 cm -2 (for N max > 10 19 cm -2 ) 
due to the steep column density profiles at these column densities 
resulting from ionization. The highest curve has N m i n equal to 
that in Rosenberg & Schneider (2003). The range of log(NHi) in 
this curve shown here is fairly narrow, thus showing that some 
variation is possible in N max for observed galaxies. 



sumed for the objects simulated in Figs. [3] through |H| In 
Fig. |3 the column densities are averaged, so that the mass 
contained within a radius where Nhi = N m in is divided by 
the area within this radius. Thus the right side of the plot 
looks similar to Fig. Q but the left side shows where the 
averaged column densities become lower when most of the 
mass in a galaxy or HI cloud is close to the ionization edge. 

Samples of 200 objects, having uniformly distributed 
values of log N ma x, are simulated in Figs. |3]and|31 We sim- 
ulate only galaxies having h = 2 kpc in Fig. |SJ as the 
curve will otherwise become widened in the steeper parts 
(N max & 10 20 cm -2 ) due to variations in disc scale lengths. 
The top curve shows objects observed to the same limits as 
in Rosenberg & Schneider (2003). It can be seen that less 
than one order of magnitude of variation in (Nhi) corre- 
sponds to more than two orders of magnitude in possible 
N m ax values. 

In Fig. [I] we plot inferred column densities (N HI ), also 
versus central column density N ma x, for simulated objects 
which are exposed to an ionizing background. Inferred col- 
umn densities are found by dividing the total HI mass of 
a galaxy or cloud by an area within some estimated HI ra- 
dius. One might use a radius where (N m in ~ 10 20 cm -2 ) 
which would be typical of the galaxies discussed in Salpeter 
and Hoffman (1996), but it would be necessary to use a ra- 
dius corresponding to a smaller N m in to find a smaller N HI . 
Thus we assume here (and for further calculations of N HI ) 
that Nmin = 10 20 cm' 2 or that logiV m i n = logJV mair - 0.5 
if log Nmin < 20.5. In Fig. 2] we show several scale length 
values between 0.2 and 4 kpc. The upper parts of the 
curves (having log N ma .x > 20.5) are similar to the central 
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Figure 4. Inferred neutral column densities are plotted versus 
central neutral column densities (N m ax) for simulated galaxies 
which are exposed to an ionizing background. Column densities 
are found here by dividing total HI masses by the area within the 
radius of N m in, where N m in = 10 20 cm -2 or that log N m in = 
log Nmax — 0.5 if log N m in < 20.5. Curves are shown for several 
disc scale length values, ranging from 0.2 to 4 kpc. Note that 
the curves are much steeper at low N max than the line shown in 
Fig. El due to ionization. 



curve in Fig. |5] although any variation with h happens only 
with an ionizing background. Here the lower curves, where 
log N m in = log Nmax — 0.5 are steepened due to ionization, 
when compared to the line with slope 0.745 without ioniza- 
tion. 

It can be seen in Figs.[j|]and[4]that the averaged and in- 
ferred column densities both fall off quickly when N max goes 
below some value ^ 10 20 cm -2 , though the N^j or (Nhi) 
value where this happens depends upon the assumed value 
of N m i n ■ The steepening of the dotted curve in Fig. [3] is a 
result of galaxy discs having ionization edges at a few x 10 19 
cm' 2 . For the inferred column densities shown in Fig.|J|the 
steepening happens in part because of the radius withinin 
which the column density is averaged, as it is for the dot- 
ted line shown in Fig. [5] However, the line in Fig. 2 has a 
slope of 0.745, whereas in Fig|l] Nfu changes by about four 
orders of magnitude when N m ax changes by two orders of 
magnitude, which is as steep as a line with a slope of 2, again 
resulting from ionization. For objects having low N m ax, the 
area of the disc which has a high column density becomes 
small, so that the HI fluxes for these objects are also likely 
to be small. 

While Nffj values may difficult to estimate, we ulti- 
mately want to know which central column density Nmax 
values can be detected in a survey having some limiting flux. 
In Figs.[!)]through[H] we plot peak fluxes for simulated galax- 
ies, which can be compared with the HIDEEP limiting value 
of 18 mjy per velocity channel, where the velocity resolution 
is 18 km/s and the channel separation is 13.2 km/s. corre- 
sponding to the vertical line in each Figure. Samples of 5000 
objects, having total gas masses between 10 6 and 10 11 M 
which obey a total gas mass function with a slope of -1.3, 
are simulated for each case. 

In Fig. central column density values are plotted ver- 



sus peak flux values which would be produced by completely 
unionized gas. More objects have lower fluxes than higher 
ones at any given Nmax, as there are more objects which 
are at larger distances, and more having lower masses. In 
Case A (black circles), where the simulated masses are cor- 
related with Nmax, and thus low Nmax objects tend to have 
smaller scale lengths, an observer would not expect to de- 
tect many objects having Nmax Zz, 10 20 cm -2 , even without 
considering the effects of ionization. Yet if the objects with 
low Nmax are as massive as those having high Nmax (Case 
B, grey circles), then numerous objects having low Nmax 
would still be seen above a reasonable flux limit. The fluxes 
for objects having lower N m ax values would only be lower 
if a substantial fraction of the mass is outside the beam, as 
seen for Nmax 10 19 cm -2 , or if some of the gas is ionized, 
as seen in Fig.|S| In Fig.|S| values of N ma x are plotted versus 
peak flux for each object, now including the effects of ioniza- 
tion. At a limiting flux above that in the HIDEEP survey, 
few objects could be detected having N m ax > 10 20 cm~ 2 
(although the detectable inferred column densities could be 
lower than this value, as seen in Fig. [HJ. It can be seen 
that an observer would need to look at least two orders of 
magnitude fainter than the HIDEEP limit to detect ionized 
objects having N ma x Ss 10 20 cm -2 in Case B, or possibly 
more if low column density objects are less massive. The 
fluxes shown, for example, for ionized objects here, may be 
lower limits if the gas far from galaxies is clumpy. 

In Figs. [7] and |H| we attempt to calculate inferred col- 
umn densities, which are plotted again versus peak flux, and 
can be compared with the N^jj values for the HIDEEP sur- 
vey. Again we find N^ T values within radii where log N m in = 
20 or log Nmin = logiVmai — 0.5 if log Nmax < 20.5, which 
makes sense, as all the objects detected in HIDEEP were 
found to have possible optical counterparts, and are thus 
assumed to have measurable HI radii. In Fig. |7] we show 
the total inferred column densities (N° ot ), or the neutral gas 
that would be seen if there were no ionizing background. A 
flat distribution of Nmax values has been assumed all along 
for Case B, while even more objects with low N ma x are sim- 
ulated in Case A, yet the majority of the points which are 
above the limiting flux are seen at N° ot ~ 10 20 ' 6 cm -2 . Thus 
the Nxj values found in HIDEEP, where the distribution is 
peaked at fo 20 65 cm -2 , should have a strongly peaked dis- 
tribution simply as a result of averaging exponential profiles 
over a radius where N m in = 10 20 cm -2 . Yet numerous lower 
N°ot would still be seen, especially for Case B. 

The Nx! values, as seen when the gas is exposed to an 
ionizing background, are shown in Fig. |H| Here only a few 
galaxies or clouds having Ngj < 10 20 cm -2 are seen above 
the HIDEEP flux limit. In either of the cases, it should be 
possible to detect more objects having low N^ T in a survey 
which is just slightly deeper than HIDEEP, if the radii of the 
objects can be estimated accurately enough. For Case A the 
objects shown appear to have a similar relationship between 
peak flux and column density as in the unionized plot, but 
the main difference is that more of the objects have inferred 
column densities that are below the minimum value shown 
on the plot. Here the distribution of detectable N^j values 
appears to be somewhat less strongly peaked compared to 
the unionized cases, yet the binned points, which are simu- 
lated like those above the HIDEEP flux limit in Fig. [HJ do 
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Figure 5. Without an ionizing background, central column den- 
sities N ma x are plotted versus peak fluxes for simulated galaxies 
having purely exponential profiles, both for a sample of objects 
where the HI masses are correlated with N m ax (Case A, black 
circles) and for a sample where the masses are independent of 
central value N max (Case B, grey circles). A line is drawn at the 
limiting flux of the HIDEEP survey. If the lower column density 
galaxies were unionized, as shown here, galaxies having a wider 
range of central column densities would be detectable. 
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Figure 7. Inferred column density values (N° ot ), without an 
ionizing background are plotted versus peak fluxes for simu- 
lated galaxies, again for the cases described in Fig. [5] assum- 
ing radii where log N m in = 20 or logN min = log N m ax - 0.5 if 
log N ma x < 20.5. The majority of the points, which are above the 
HIDEEP limiting flux, are seen at 10 20 ' 6 cm -2 . Some lower N° ot 
values could be detected, in this hypothetical unionized case, if 
the radii of these objects were correctly estimated. 
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Figure 6. Central column density (N m ax) values are plotted 
versus peak fluxes for simulated galaxies, for the cases shown in 
Fig. |5] but now exposed to an ionizing background. Few points 
having N m ax below 10 20 cm -2 can be detected below the limiting 
flux of the HIDEEP survey. 
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Figure 8. Inferred column density (Nfrr) values are plotted ver- 
sus peak fluxes for simulated galaxies, for the two cases described 
in Fig. [5] Here the gas has been exposed to an ionizing back- 
ground, and N m i n values are assumed as in Fig. |7| Above the 
limiting HIDEEP flux, only a few points might be detected hav- 
ing Nmax & 10 20 cm" 2 . 



not look very different from the HIDEEP distribution, as 
shown in Figs. Inland I1UI 

In Figs , 19*1 and 1 1 01 we show the inferred column density 
distributions for Cases A and B respectively, for samples 
of 1000 objects having peak fluxes > 18 mjy and velocity 
widths > 40 km/s. Also shown is the Gaussian curve fit- 
ted to the HIDEEP galaxies, having a mean of 20.65 and a 
scatter of 0.38, as in Minchin et al. (2003) and the binned 
data from Minchin et al. (2003). Neither histogram is very 
different from the Gaussian or data set, assuming some mea- 
surement uncertainties. Case A appears to be somewhat 



strongly peaked here, although the observational uncertain- 
ties, mostly in measuring the area of the galaxies, are not yet 
shown here. Case B appears to be more broad, and similar 
to the plotted Gaussian. The peak feature at 10 20 ' 3 cm -2 
is also seen in the Minchin et al. (2003) data. This feature 
is a result of many galaxies with high N ma x having their 
column densities averaged over similar Nmin, and confirms 
that similar galaxy radii are being used to find the inferred 
column densities for the majority of the galaxies for the ob- 
servations and for the simulations done here. The Case A 
model looks more realistic at the high column density end 
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Figure 9. Distributions are shown for simulated inferred column 
density (iVjj-r) values detected above the limiting flux of HIDEEP, 
for Case A, where the gas mass for each galaxy is closely related 
to Nmax (dashed line). Also shown is the Gaussian curve fitted 
to the HIDEEP sample (solid line) and the binned HIDEEP data 
(dotted line). 



Figure 10. Distributions are shown for simulated iVj^, values 
detected above the limiting flux of HIDEEP, for Case B, where 
the total gas mass is independent of Nmax (dashed line). Also 
shown is the Gaussian curve fitted to the HIDEEP sample (solid 
line) and the binned HIDEEP data (dotted line). 



of the distribution, although we do not model high column 
density galaxies carefully given that more of their gas may 
be converted into stars or ionized by these stars. Neither dis- 
tribution is very different from the Gaussian curve, although 
the most realistic scenario is likely to have an intermediate 
behaviour between Cases A and B. There is likely to be some 
relationship between the gas masses and central column den- 
sities of galaxies, but the scatter may be larger than we have 
assumed in Case A. For example, we increase the amount 
of scatter for N max from 0.5 to 1.5 orders of magnitude in 
Case C. 

In Case B (Fig. [TUJl there are a few objects seen which 
have sufficiently high fluxes to be detected, but column den- 
sities below anything detected by HIDEEP. (These points 
are not seen for Case A, only because sufficiently low column 
densities are assumed to arise only in very low mass galax- 
ies which are all below the HIDEEP flux limit.) Such points 
may have simply not yet been detected by HIDEEP due to 
the small number of objects detected, compared to the 1000 
objects simulated here. Also, these objects may have small, 
and thus uncertain radii. Very few HI clouds are thought to 
have no optical counterparts (Davies et al. 2004), and there 
are no clouds detected which have no optical counterparts 
in HIDEEP. This could happen because the gas is actu- 
ally more clumpy than we have modelled here, so that star 
formation occurs in small regions within these clouds. The 
isolated HI cloud detected by Giovanelli & Haynes (1989), 
for example, was later found by many observers to contain 
a small dwarf galaxy. The radii for such objects are thus 
likely to be underestimated, so that their inferred column 
densities will be higher than what we find here. 

The inferred column density distributions for the union- 
ized cases are surprisingly indistinguishable from those 
shown in Figs. 151 and 1101 when assuming that galaxies hav- 
ing velocity widths < 40 km/s are not detectable, although 
these objects were not removed from Figs. |^|through|S| Thus 
velocity width related selection effects could be as important 



as ionization in determining the shape of the distribution of 
inferred column densities. Objects which could be ionized to 
the point of being undetectable have low masses and thus 
low velocity widths in Case A. In Case B the distribution 
of Nffj is the same for detected and undetected objects. 
However the neutral and ionized column density distribu- 
tions are not the same in the intermediate Case C, where 
\ogN max — 21.7+1.0 log(M tot — 10. 0)±1. 5, as some galaxies 
having detectable velocity widths are affected by ionization, 
yet there is some variation with galaxy mass. In Case C the 
iV t ° ot distribution is very strongly peaked, yet the N^j dis- 
tribution is somewhat broad, and similar to that seen for 
Case B. 

The biggest source of uncertainty in finding inferred col- 
umn densities is in measuring the galaxy radii. When an un- 
certainty of 40% in the radii is included, as shown in Fig. II II 
Case A is still too strongly peaked, while Cases B and C are 
more similar to each other and to the Gaussian curve fit- 
ted to the HIDEEP data. Fluctuations in the ionizing back- 
ground radiation are also likely to broaden the distributions 
somewhat, but only by increasing the number of galaxies 
with high Nfjj as there are few locations in space (Linder 
et al. 2003) where the ionizing background is as low as the 
value assumed here (Haardt & Madau 1996). 

A summary of the simulations is given in Table[T| where 
the mean and scatter values for the Ngj distributions are 
listed. Case C is intended to be intermediate between the 
Cases A and B shown previously, as there is some correlation 
between Nmax and the gas masses, but a larger scatter so 
that logN max = 21.7+ 1.0 log(M tot - 10.0) ± 1.5. 

In an attempt to put some constraints on the properties 
of the objects detected in the HIDEEP survey, some further 
variations were made of the model parameters. The slope in 
the log Nmax-log Mtot relationship was increased from 1.0 
to 1.2, which resulted in a strongly peaked inferred column 
density distribution similar to that for Case A. In further 
simulations we assume that log N ma x = 21.7+1.0 log(Mt t — 
10.0) ± 1.5 as in Case C. Flattening the central parts of the 
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Table 1. Inferred Column Densities and Lyman Limit Absorber 
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Figure 11. Distributions of simulated N^j (dashed lines) are 
shown with uncertainty in the measured galaxy radii included 
for Cases A (left), B (centre), and C (right). Also shown is the 
Gaussian curve fitted to the HIDEEP sample in each frame (solid 
line). 



galaxy column density profiles (to a value between 10 20 ' 5 
and 10 21 ' 5 cm" 2 ) also gave rise to a strongly peaked distri- 
bution as in Case A. We also simulated a stronger ionizing 
background, based upon the z < 1 measurements of Scott 
et al. (2002), and did a simulation using a steeper slope 
of —1.5 for the HI mass function, as reported by Rosen- 
berg & Schneider (2002). In both of these cases the in- 
ferred column distribution is similar to that in Case C. Fur- 
thermore we varied the characteristic Nmax value, so that 
log AW = 22.2 + 1.01og(M to i - 10.0) ± 1.5 while using 
flattened central column density profiles in order to allow 
for higher extrapolated N max values as suggested in Bowen, 
Blades, & Pettini (1996). In this case the Nfjj distribution 
is somewhat strongly peaked. A further simulation was done 
where the exponential column density profiles were flattened 
to a power law with a slope of —4 beyond four HI scale 
lengths, in which case we see a slightly excessive number 
of high column density galaxies. Combining the slow fall-off 
in the outer parts with a reduced central column densities 
would likely give rise to a more realistic number of high col- 
umn density galaxies. For any of the parameter variations, 
the distribution is not very different from the Gaussian curve 
fitted to the HIDEEP points, given the uncertainty in the 
actual relationship between gas masses and central column 
densities. 



5 LYMAN LIMIT ABSORBER COUNTS 

Lyman limit absorber counts, arising from quasar lines of 
sight through gas having 10 17 ' 2 < N H i < 10 20 3 cm" 2 pro- 
vide further constraints on the numbers and properties of 
undetected objects containing low column density gas. Esti- 
mates have been made for the number of Lyman limit sys- 
tems arising in optically observed galaxies, but this gener- 
ally involves assuming a cross section for absorption around 
a galaxy as a function of its optical luminosity (Linder et 
al. 2003; Steidel 1995; Bergeron & Boisse 1991). However 
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Mean inferred column densities, scatter for the inferred column 
density distributions, and Lyman limit absorber counts are shown 
for the simulations as summarized here: 

Case A, where gas masses are related to Nmax such that 

log AW = 21.7 + 1.0 log(A/ tot - 10.0) ±0.5; 

Case B, where gas masses are independent of Nmax', 

Case C, where logAWc = 21.7 + 1.0 log(Af tot - 10.0) ± 1.5; 

Cases Au, Bu, and Cu, versions including uncertainties in galaxy 

radii. 



the relationship between optical and HI properties of galax- 
ies may not be well enough understood to make such es- 
timates. Here we can estimate the number of Lyman limit 
systems arising from galaxies that obey an observed HI mass 
function instead. 

We estimate the number of Lyman limit systems arising 
in each scenario by putting random lines of sight through the 
sphere in which the galaxies are simulated and finding the 
column density where a line of sight intersects a galaxy. We 
simulate 20,000 galaxies in an eighth of a sphere having a 
radius of 108 Mpc, and use 10,000 lines of sight for each case. 
The number of Lyman limit systems is then calculated by 
correcting to a number density of simulated objects which 
gives rise to an HI mass function having a normalization 
consistent with Zwaan et al. (1997). The minimum galaxy 
mass of 1O 7 M0 is used for simulated galaxies, as the lowest 
mass objects are not likely to make a substantial contri- 
bution to Lyman limit (or lower column density) absorber 
counts (Linder 1998), although estimates from HI studies 
suggest that slightly more massive objects do contribute to 
Lyman limit absorption (Ryan- Weber et al. 2003). Values 
for the number of Lyman limit absorbers per unit redshift 
along a line of sight, (dN/dz)o,LL, are shown for the three 
main simulations in Table 

The number of Lyman limit absorbers has been mea- 
sured at redshifts ^ 0.36, and the evolution is seen to 
be about flat or slightly decreasing down to redshift zero. 
The lowest redshift values available are within the range 
of (dN/dz) LL ~ 0.2 to 1.3 (Lanzetta, Wolfe, & Turnshek 
1995b; Storrie-Lombardi et al. 1994; Stengler-Larrea et al. 
1995). Most of the simulations appear to be consistent with 
these observations, although Case A gives rise to too many 
absorbers. The simulation where the ionizing background in- 
tensity was increased gives rise to too few (0.08) absorbers 
per unit redshift, but the ionizing intensity used is proba- 
bly more relevant at z ~ 1, where there are actually more 
absorbers because less cosmological expansion has occurred. 

Case A (and some other simulations giving more 
strongly peaked N%j distributions) appear to give rise to a 
somewhat excessive numbers of Lyman limit systems. How- 
ever, the same problem seems to arise when estimating the 
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number of Lyman limit systems around galaxies whose op- 
tical luminosity function is known. It has long been thought 
that luminous galaxies have sufficient cross sections to ex- 
plain the Lyman limit absorber counts fully, yet it is not 
known why dwarf and LSB galaxies would not also make 
some contribution, especially now that such faint objects 
are often found to give rise to damped Lya absorption. It 
is possible, for example, that feedback processes change the 
column density profiles in the outer parts of some galaxies 
(McLin, Giroux, & Stocke 1998). While we do not rule out 
cases simply because they give rise to somewhat excessive 
Lyman limit absorber counts, Case A also has an excessively 
peaked distribution of inferred column densities, which is 
not improved when taking uncertainties in measuring the 
galaxy radii into consideration. Thus allowing for a wider 
range of galaxy column density profiles could make Lyman 
limit absorber counts more consistent with what we expect 
from observed galaxies. 



6 CONCLUSIONS 

Most galaxies have inferred column densities around 
10 20,6±0 ' 3 cm~ 2 because inferred column densities are found 
by averaging column density profiles, which are exponential 
or similar, over a radius where the minimum column den- 
sity is ~ 10 20 cm -2 . Ionization plays some role in making 
lower column density objects undetectable, including those 
without substantial optical counterparts. However inferred 
column density distributions tell us little about the distri- 
bution of central column densities in galaxies and clouds. 

Ionization by the background of ultraviolet photons will 
strongly affect the amount of neutral gas remaining, and 
thus the HI flux detected, in objects having low hydrogen 
column densities, if such objects having sizes comparable to 
galaxies exist. Typical HI fluxes are reduced, as a result of 
ionization, by a factor of ~ 100 for galaxies having peak col- 
umn densities N max ~ 10 19 ' 5 cm~ 2 compared to those with 
N m ax ~ 10 20 cm~ 2 , even if the lower column density galax- 
ies are extended in size and just as massive as the higher 
column density galaxies. 

We do not always know the central column densities 
of the faintest HI sources, but the detected inferred col- 
umn densities are also likely to be above ~ 10 20 cm~ 2 
for most observable galaxies. Inferred column densities are 
rather weakly related to central column densities for objects 
having exponential profiles. Furthermore, since HI profiles 
tend to be mapped out to limiting column densities ~ 10 20 
cm~ 2 , it may be difficult to estimate the radii, and thus 
the inferred column densities in a consistent manner for ob- 
jects having lower N max values. For example, if the radii 
are underestimated, which might be more likely to happen 
for an extended, diffuse galaxy, the inferred column den- 
sity could be overestimated. Other selection effects, such as 
those against objects having low velocity widths, may also 
be important in understanding the observed distribution of 
inferred HI column densities. 

The observed distribution of inferred HI column den- 
sities, as seen by Minchin et al. (2003), can easily be sim- 
ulated assuming possible populations of galaxies having a 
wide range of size and central column density distributions, 
and the simulated distributions arc similar to the HIDEEP 



distribution for a wide range of model parameters. (Thus 
the 'Frozen Disc' hypothesis of Minchin et al. 2003 seems 
to be unnecessary in explaining these observations.) How- 
ever, we are thus given little constraint on the properties of 
gas rich objects which have so far escaped detection in the 
deepest HI surveys. Given the effects of ionization, we are 
unable to rule out the existence of undetected populations 
of very faint dwarf galaxies or giant gas clouds, as long as 
they have low central column densities. Such objects could 
make some contribution to Lya absorption, although a more 
reasonable number of Lyman limit systems arises if galaxies 
have a wide, rather than narrow, range of central column 
densities. 

The ionizing background radiation is more intense at 
redshifts around 1 or 2 than at redshift zero (Haardt & 
Madau 1996), and therefore some of the apparently younger 
galaxies, such as LSB galaxies, may have been ionized at 
these redshifts if they have lower central column densities 
(de Blok et al. 1996), thus slowing their evolution. Ionization 
may have also affected the formation of dwarf galaxies in cer- 
tain environments at high redshifts (Efstathiou 1992; Tully 
et al. 2002), as less dense environments are more likely to be 
optically thin to ionizing radiation when the dwarf galaxies 
formed. Thus dwarf galaxies may have formed more easily 
in rich clusters such as Virgo (Sabatini et al. 2003) and For- 
nax (Kambas et al. 2000) than in more diffuse clusters such 
as Ursa Major (Trentham & Tully 2002) and other envi- 
ronments (Roberts et al. 2003) . Understanding the role that 
ionization plays is thus important in testing Cold Dark Mat- 
ter scenarios and other theories related to galaxy formation. 
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